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Abstract 

We use spin-density-functional theory within an envelope function approach to calculate elec¬ 
tronic states in a GaAs/InAs core-shell nanowire pierced by an axial magnetic field. Our fully 
3D quantum modeling includes explicitly the description of the realistic cross-section and compo¬ 
sition of the sample, and the electrostatic field induced by external gates in two different device 
geometries, gate-all-around and back-gate. At low magnetic fields, we investigate Aharonov-Bohm 
oscillations and signatures therein of the discrete symmetry of the electronic system, and we criti¬ 
cally analyze recent magnetoconductance observations. At high magnetic fields we find that several 
charge and spin transitions occur. We discuss the origin of these transitions in terms of different 
localization and Coulomb regimes and predict their signatures in magnetoconductance experiments. 
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I. INTRODUCTION 


Gated semiconductor nanowire (NW) devices represent flexible test beds to study trans¬ 
port phenomena in the quasi-ID quantum regime. In this context, InAs-based NWs offer 
privileged properties derived, for instance, from the light InAs electron effective mass, which 
allows the experimental observation of the subband spectrum quantization even in NWs of 
relatively large section , 1 3 or from its large spin-orbit, interaction and Lande factor ] 45 which 
boost their prospectives in spintronics® even at relatively high temperature ) 7 Furthermore, 
in this narrow gap material the Fermi energy, Ep, is pinned by surface states above the 
conduction band edge® leading to an accumulation of electrons at the NW surface and 
facilitating the fabrication of ohmic contacts.®^ 

The resulting tubular shape of the conducting channel points toward interesting quan¬ 
tum phenomena under external magnetic fields ! 10 In particular, an axial field may lead to 
Aharonov-Bohm (AB) field-periodic modulation of the electron energy spectrum !- 11 And, if the 
phase coherent length exceeds the perimeter of the NW, the observation of magnetoconduc¬ 
tance oscillations . 12 13 Indeed, several observations of AB-like oscillations in magnetotrans¬ 
port experiments performed on radial heterostructures have been reported . 14 18 Recently, 
Girl et a /. 16 have observed flux-periodic magnetoconductance oscillations in GaAs/InAs 
core-shell NWs. The oscillations persisted at different density regimes, modulated by a 
back-gate, exhibiting phase shifts as the back-gate voltage was gradually increased. A field- 
periodic magnetoconductance has also been observed in the same system with superconduc¬ 
tor contact ^ 17 and, after removal of the GaAs core, in a hollow InAs shell . 18 

The single-crystal NW-based heterostructures investigated in these experiments have a 
prismatic hexagonal cross-section. However, the experimental observations were analyzed in 
terms of simplified cylindrical electronic systems, and the potential induced by the back-gate 
voltage, which also removes the cylindrical symmetry, was neglected. Likewise, theoretical 
calculations dealing with radial electronic systems with an axial magnetic field usually as¬ 
sume a cylindrical symmetry . 5 12 13 19 20 Ferrari et al . 10 investigated the effect of an axial 
magnetic field in prismatic systems, but the single-particle model adopted did not allow for 
a direct comparison with experiments. 

Such approximations are particularly severe in radial heterostructures, where coupling 
between the discrete (hexagonal in InAs or GaAs) symmetry and many-electron interactions 
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leads to strongly inhomogeneously distributed electron gas and, in turn, to the coexistence of 
ID and 2D channels at the corners and facets of the hexagonal hetero-interfaces.® 23 Strong 
anisotropy-induced effects are predicted in this case, such as negative magneto-resistance in a 
transverse magnetic field 24 and symmetry-induced cancellation of the AB effect in hexagonal 
quantum rings 25 The inhomogeneous electron gas localization was crucially exposed in the 
recent observation of intra- and inter-band excitations. 23 26 

In this paper we study the electronic states and magnetoconductance in GaAs/InAs 
core-shell NWs with an axial magnetic field within a spin-density-functional theory (SDFT) 
approach. Our fully 3D modeling explicitly includes the description of the quantum states 
within an envelope function approach with realistic cross-section and composition of the 
sample, and includes the electrostatic field induced by external gates in two different de¬ 
vice geometries, gate-all-around and back-gate. At low magnetic fields, we investigate the 
nature of the magnetoconductance oscillations, as measured in Ref. [161 predicting specific 
signatures of the discrete symmetry of the electronic system in the AB magnetoconductance 
oscillations, and justifying the observation of AB oscillations despite the broken symmetry 
induced by the back-gate voltage. At high magnetic fields we found that several charge 
and spin transitions occur. We discuss the origin of these transitions in terms of different 
magnetic field induced localization and Coulomb regimes and predict their signatures in 
magnetoconductance experiments. 


II. THEORETICAL MODEL 

Within a parabolic single-band envelope-function description, the effective Kohn-Sham 
Hamiltonian under an external magnetic field reads 

H = l (P — e A(R)) (P — e A(R)) + V conJ (R) + V/(R) + V H (R) + V} C (R). (1) 

Here, R = ( x,y,z ), P is the momentum operator, A(R) is the vector potential, e is the 
elementary charge, and m*( R) is the material-dependent electron effective mass. V con f(TV) is 
the spatial confinement potential induced by the heterostructure and Vfl-(R) is the Hartree 
potential energy. The Zeeman energy (R) and the exchange-correlation potential V^ci R) 
depend on the the spin index a =j\ i of the electrons. 
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We consider an infinitely long NW extending along the z direction. To describe an 
axial magnetic field we adopt the symmetric gauge A(R) = B/2(—y,x,0) (see Fig. 1 for 
axis definition). The axial field does not break the spatial invariance along the z-axis. 
Therefore, the single-particle eigenfunctions of Q can be written as T nj fc j(T (R) = e* fcz 0 niCT (r), 
with r = (x,y), n the principal quantum number and k the wavenumber along direction z. 
Substituting ^„,A:,ct(R) and A(R) in (Jl]) we obtain the spin-dependent Kohn-Sham equation 


h 2 v 1 y 

2 r m*(r) r ' 2m*(r) z ' 8m*(r) 


eB f e 2 B 2 
-L z + 


x 2 + y 2 ) + v conf (r) 


+4( r ) +un(r) +< c (r) 


0n,<r(l’) £n,k,cr 4*n,a(x')■ 


( 2 ) 


Here, t n .k. a = £ n ,a + includes the ID parabolic dispersion along the z-axis, and L z = 
~ ih {x?ij}-y&) is the azimuthal angular momentum operator. To obtain Eq. (2) it is 
necessary to assume that the ^-component of the effective mass, m* z , does not depend on r, 
i.e., on the material. This approximation is expected to have a small effectP 7 and enables to 
decouple the electron motion in the longitudinal and transverse directions. 

The confinement potential v con f( r) is set by the conduction band offsets among the dif¬ 
ferent materials that are radially modulated in the NW cross section. The Zeeman term 
is 


v z( r ) = 9 *(r)fi B Br] a , ( 3 ) 

where g*( r) is the material dependent Lande factor, ys is the Bohr magneton, and = 
+l/2(—1/2) for a =t (|). 

The Hartree potential energy, vh(f), is calculated from the electrostatic potential, 
vh( r) = —e<h(r), via the Poisson equation 

Ve:(r)V < l>(r) = — e(n( r) — n D ( r)) . (4) 

£o 

Here, n(r) = n- j-(r) + n^(r) is the total free electron charge density calculated, using the 
Khon-Sham eigenstates obtained from Eq. , as 

i _ r°° 

n v( r ) = I 2 / dk f(e nAa -E f ,T), (5) 

2 ^ _ J-oo 
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where 


f{tn,k,CT ~ Ep, T ) — 


( 6 ) 


X _|_ (A e n,k,a-E F )/k B T 1 

is the Fermi occupation, with Ep, T and ftp being, respectively, the Fermi energy, temper¬ 
ature, and Boltzmann constant. In Eq. Q np>( r) is the density of static donors and e(r) is 
the material dependent static dielectric constant. 

The exchange and correlation potential, tE,(r), in the local-spin-density approximation 
(LSDA) is given by the functional derivative 


= 6e xc K r),C(r)) 
xc[ } Sn a (r) 

where e xc (n(r), C( r )) is the exchange and correlation energy density and 

^ \ _ n t( r ) - ™t( r ) 


(7) 

( 8 ) 


is the local spin polarization. In the present paper, we use the correlation functional proposed 
by Per dew and Wang. 28 

From the solutions of the Kohn-Sham equations we also obtain the total free energy per 
unit length froml 29 


E = — 


1 

2~7T 


[ dk fn,k,a — \.e f drv H (r)n(r) - ^ / dr 

n,tr Z J a J 


V xc( T ) n <r( r ) 


kn T 


T j dr £ xc (?7.(f), (((l - )) T ^ ^ j dk [fn,k,cr in fn,k,cr T (1 fn,k,rr') ln(l /n,fc,cr)] ■ 


71,(7 


(9) 


Here, the second term on the right hand side is the Hartree energy per unit length with the 
sign inverted, the fourth term is the exchange and correlation energy per unit length, and 
the last term is an entropy functional, where f n ,k,a = f(e n ,k,a ~ Ep, T ). 

Equations (§-@i are solved iteratively until self-consistency is reached, which we con¬ 
sider to occur when two convergence criteria are simultaneously fulhlled in two consecutive 
iterations: first, the relative variation of the charge density is lower than ICE 4 at any point 
of the discretization domain, and second, the relative variation in total free energy per unit 
length (Eq. [9]) is lower than ICE 8 . 
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Equations (|2]) and Q are numerically integrated in a real space hexagonal domain. We 
use the same symmetry preserving triangular grid with ~ 1.14 points/nm 2 for both for¬ 
mulas and integrate Eqs. [2] and [4] with the methods of finite elements and finite volumes, 
respectively. Dirichlet boundary conditions are assumed in both cases, generally forcing 
the solutions to vanish at the boundaries. To simulate the effect of a gate-all-around (see 
Rg.0 a)), the electrostatic potential in the Poisson equation is forced to take the gate volt¬ 
age V g at the domain boundaries. For a back-gate, we assume that the hexagonal domain 
is sandwiched by two flat infinite electrodes (see Fig. [l](b)) and the electrostatic potential 
is set at the gate voltage V g at the bottom facet and zero at the top one. Accordingly, at 
the lateral boundaries the electrostatic potential is set to F ■ ds(r), with F and de(r) being, 
respectively, the electric field in the capacitor and the vertical distance from the boundary 
point to the bottom electrode (see Fig. [jjb)). 

Finally, we also calculate the spin-projected free charge density per unit length 


n 


a 



( 10 ) 


and the spin-projected ballistic conductance of the NW by means of the linear-response 
Landauer formula 




df{E — E Fl T) 
dE 


dE , 


( 11 ) 


where the integral is performed along each energy spin-subband B n)t7 . Note that the inte¬ 
grand gives a significant contribution only in the energy region close to the crossings of the 
subbband with the Ferrny energy Ep. 


III. NUMERICAL RESULTS 

We consider a GaAs/InAs core-shell hexagonal NW as the one measured in Ref. [16] and 
outlined in Fig. [lj It is composed of a GaAs NW core with a minimal diameter of 100 
nm and an InAs shell with a thickness of 25 nm. In addition, we include in the device an 
external 30-nm-thick layer of Si02 intended to simulate the insulating layer that separates 
the conducting channel from a back-gate in the experiment A 6 The GaAs core is doped with 
a homogeneous density of donors rip, = 5 x 10 15 cm -3 . The material parameters used in 
the calculations are listed in Tab. |T| where the conduction band edge, Ecb, is obtained with 
the so-called 40:60 rule 30 31 from the band gap. Calculations have been conducted assuming 
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FIG. 1. Schematics of a core-shell NW in the (a) gate-all-around and (b) back-gate device config¬ 
urations. 



GaAs 

InAs 

Si0 2 

sk 

m 

0.067 

0.028 

0.41 

£ 

13.18 

15.5 

3.9 

9* 

-0.484 

-14.3 

2.0 

Ecb i 

(eV) 0.858 

0.252 

5.4 


TABLE I. Material parameters used in the simulations; electron effective mass (m*), dielectric 
constant (e), effective Lande factor (g*) and conduction band edge ( Ecb )• 


a Fermi energy placed 75 meV above the InAs conduction band edge (as in Ref. H6j), a 
temperature of 1.8 K, and the InAs effective mass as the constant mass factor (m*) arising 
in the parabolic dispersion of the ID subbands (see Eq. i§). 
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FIG. 2. (a) Magnetic spin-subbands (MSS) in the low field regime. Red and blue dots indicate f- 
MSSs and "f-MSSs, respectively. The horizontal black line is set at Ep. (b) Self-consistent electron 
density distribution, n(r), for the InAs/GaAs NW at B = 0. (c) Spin-projected magnetoconduc¬ 
tances Gy (red) and G ^ (blue), and total magnetoconductance (black). 


A. Low magnetic field regime: magnetoconductance oscillations 

In Fig.[2]we show the ground state properties and magnetoconductance of the investigated 
core-shell NW at V g — 0. The density distribution of conduction band electrons (Fig. [2j)b) ) 
shows that charge is exclusively accumulated in the InAs shell and preferentially localized 
at the corners of the hexagonal section. As reported for several core-(multi)shell hexagonal 
NWsJ 2H2M2J such distribution is favored by Coulomb interactions which tend to increase 
the inter-electron distance. In Fig. [2j)a) we show the energies of the spin-subbands edges at 
different magnetic fields, hereafter referred as magnetic spin-subbands (MSS), with spin up 
(t-MSSs) and spin down (J,-MSSs). Due to the hexagonal symmetry of the self-consistent 
potential, the low-energy spectrum is at low fields formed of groups of 12 MSSs arising from 
the 6 irreducible representations of the C6 symmetry group. Each of these groups is further 
spin-split by the strong Zeeman effect forming two bunches, one of f-MSS and the other of 
4,-MSSs, of 6 braided MSSs. 

Within each group the MSSs oscillate due to the AB effect, developing crossings with 
MSSs of their same group, that have different symmetry and/or different spin, and anticross- 
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ings with MSSs of neighbouring groups with same symmetry and spin. The oscillation period 
is ~ 0.32 T. Since the calculated expectation value of the radial position, p = \Jx 2 + y 2 , 
of the electron system is 66.36 nm, this periodicity fairly corresponds to the periodicity of 
~ 0.30 T of the corresponding circular system. 

In Fig. gc ) we show the spin-projected magnetoconductances G a (B) and the total mag¬ 
netoconductance G(B ) = G^(B) + Gi(B). Even though both G a (B ) exhibit regular flux- 
periodic oscillations, G(B ) only does so at very low fields. After the second oscillation cycle 
the G(B) periodicity is suppressed by the Zeeman effect which breaks the periodicity of the 
MSS spectrm. 12 13 Apart from this, G{B) does not differ qualitatively from that of an elec¬ 
tron system in a cylindrical tube! 12 13 16 Indeed, in the present case, Ep lies within one group 
of braided MSSs, and the spectrum around Ep is similar to that of a cylindrical system. 
However, in an experiment Ep can be tuned by means of external gates. Therefore, we next 
study the system at different Fermi levels Ep or applied gate voltages V g . 

In Fig. [3] we show the effect of a gate-all-around voltage. This geometry tunes the 
position of the MSSs with respect to Ep , modulating the total density in the system while 
preserving the hexagonal symmetry. As shown in Fig. ga), the oscillatory behavior of the 
magnetoconductance due to the AB effect is absent at certain voltages. For instance, at 
V g = 80 mV the magnetoconductance is flat. This is due to the positioning of Ep in the 
energy gap between the second and third group of MSSs, as shown in Fig. gc). Since 
Ep does not cross any MSS, the number of conducting channels is constant. Comparing 
Figs. gb),(c), which correspond to V g =-60 and 80 mV, respectively, we also observe that 
the gate voltage affects both the width of the MSSs groups and the gaps between them. In 
fact, V g affects the total electron density and, hence, electron localization. As shown in the 
insets of Figs. gb),(c), a V g > 0 favors localization in the corners of the InAs shell, due to 
the larger electron-electron interaction. This, in turn, reduces the tunneling among states 
at the corners and, hence, the splittings within bunches of MSSs, while it increases the gaps 
between consecutive bunches^ 0 Note that, since the latter gaps are a direct consequence of 
the discrete symmetry of the system, flat magnetoconductance is a direct signature of the 
hexagonal symmetry, which is more likely to be observed at positive gate voltages. 

Observation of flat magnetoconductance when sweeping V g has not been reported in 
the transport measurements performed hitherto on hexagonal NWs under axial magnetic 
fields However, in these works the electron density was normally modulated by a 
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FIG. 3. (a) Total magnetoconductance at selected gate-all-around voltages, V g , as indicated by 
labels (the V g = 0 curve is the same as the black line in Fig. [2](c)) . (b) MSSs at V g = —60mV. (c) 
MSSs at V g = 80mV. Insets in (b), (c) show the corresponding n(r). 


back-gate instead of a gate-all-around. The electrostatic field generated by a back-gate 
removes the hexagonal symmetry of the electronic system and it could even destroy the 
doubly-connected topology that originates the AB effect. Therefore, one may wonder why 
flux-periodic oscillations in the magnetoconductance are observed at all. 

In order to asses this point, in Fig. [4] we show the results of simulations performed at 
different back-gate voltages. As shown in the insets of Figs. |4](b) , (c), the applied voltage 
strongly reshapes the electron density distribution in the NW. At negative (positive) V g the 
total density in the system is reduced (increased) and concentrated in the top (bottom) 
half of the InAs shell. However, whereas the doubly-connected topology which results in 
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FIG. 4. (a) Same as in Fig. [3] but for a back-gate device, (b) MSSs at V g = — 80mV (c) MSSs at 
V g = 200mV. 


AB oscillations is removed at V g < 0, it is robust for V g > 0. The origin of this difference 
can be appreciated from the corresponding MSSs (Figs. |4|b),(c)). The lowest lying MSSs 
are strongly affected by the gate, loosing the doubly-connected topology and showing an 
almost linear dispersion with the magnetic field. Higher energy MSSs, on the contrary, be¬ 
ing more delocalized over the NW section, still show doubly-connected topology. Since at 
V g < 0 only low-lying MSSs are occupied (see Fig. |4](b)) , the total electron density loses 
the doubly-connected topology and the corresponding magnetoconductance does not show 
AB oscillations. By contrast, at V g > 0 several states with the doubly-connected topology 
are occupied and the AB oscillations of the magnetoconductance persist (see Fig. 0a)). 
The latter is indeed the usual regime in magnetotransport experiment^- 16 17 where, there- 
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fore, periodic oscillations in the magnetoconductance are observed despite the symmetry 
reduction. 

B. High magnetic field regime: spin and charge transitions 




0 1 2 3 4 5 


(10 16 cnr 3 ) 


FIG. 5. (a) MSSs up to complete charge depletion. Blue (red) dots are used for f-MSSs (|- 
MSSs). The horizontal line indicates the position of Ep. Vertical dashed arrows indicate fields 
at which different spin/charge transitions occur (see text), (b)-(f) Self-consistent electron density 
distributions n(r) at selected magnetic fields. 

We next study the high magnetic held regime, up to the limit of complete electron 
depletion, which occurs at B ~ 20 T in this sample. Fig. [5] shows the MSSs and the 
self-consistent total electron density distributions at selected fields (spin-projected electron 
densities show only minor differencies and are not shown here). All simulations in this 
section are performed at V g = 0. The overall behaviour of MSSs shows that, in addition to 
the diamagnetic shift, several transitions occur at discrete fields, as we discuss below. 

The evolution of n( r) in Figs. |5](b)-(f) shows that the axial held induces a transition from 
an electron distribution localized at the corners (low held, Figs. [5j)a) ,(b)) to a distribution 
increasingly localized in the center of the facets (high held, Figs. [5j)e) ,(f)). This charge 
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reshaping is induced by the diamagnetic term (third term in the left hand side of Eq. [2]), 
which constrains the electron density to adopt distributions with lower radius as the held is 
increased, counteracted by Coulomb interactions. 

Such a corner-to-facet transition can be correlated with the evolution of the MSSs. In 
Fig. [5|a) the lowest lying bunch of 12 MSSs at B = 0 corresponds to states localized at the 
corners, whereas the second set of states are localized at the facets for orthogonality. As the 
held is increased Zeeman spin-splitting takes place and the two sets of 6 J,-MSSs approach 
in energy eventually overlapping at Bc-j.f ~ 10.2 T. At this point, the 2D electron density 
integrated along the minimal (facet-to-facet) and maximal (corner-to-corner) diameter 22 is 
nearly the same (see Fig. |5](d)) . At B > Bc_>f the 6 lowest J,-MSSs are localized at the 
facets of the inner interface whilst corner states are much higher in energy, corresponding 
to the third group of 6 J,-MSS. The same transition occurs for jvMSS, however, these states 
are already depopulated at the transition held. 

Apart from this smooth spatial localization transition, two abrupt changes of slope appear 
in the calculated MSSs. The hrst one occurs at Bp = 7.5 T and corresponds to complete 
spin polarization, as demonstrated by the spin-projected electron densities hq, n± shown in 
Fig.ga) and the corresponding spin polarization in Fig. [b](b), which marks a clear transition 
to a ferromagnetic state at B P . Note that the total density (black line in Fig. |6](a)) is reduced 
by the magnetic held with a roughly parabolic trend due to the depletion of successive, high 
energy MSSs. However, the curve shows a change of slope at (Bp). At helds right after Bp 
the rate at which the NW is depleted decreases momentarily, hq passes abruptly from being 
increased to decreased at B P , in agreement with the inversion of the J,-MMSs slope exposed 
in Fig. [5|a). 

The singular behavior of the spin polarization (Fig. |6](b)) is reminiscent of the hrst-order 
phase transition of a 2D electron gas with an in-plane magnetic held 33 34 (note that in our 
system the Seitz radius r s ~ 0.07 at zero held, which is a very weakly correlated regime), 
although it is difficult in our numerical treatment to establish whether it is a weakly hrst 
order or continuous transition. The inset in Fig. [6](b) shows the spin susceptibility, i.e., the 
magnetic held derivative of the spin polarization. This magnitude oscillates with the held 
as a consequence of the interplay between the AB effect and the Zeeman splitting which 
produce short period modulations of the spin densities. 

At helds higher than B P and Bc->f the MSSs shown in Fig. [5](a) rearrange in groups of 
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FIG. 6. (a) Total density h (black lines) and spin-projected densities n^ (blue lines) and (red 
lines) as a function of the field intensity B. Vertical dashed lines illustrate the transition fields in 
Fig. [5](a). (b) Spin polarization as a function of the magnetic field. Inset: spin susceptibility. 

six which tend to form Landau-like bands. Finally, at a larger held Bi = 16 T the spectrum 
shows an additional transition. This corresponds to complete depletion of the incipient 
second Landau-likc band. The transition is also marked by a weak but visible kink in n(B ), 
as shown in Fig. [6](a), which, as for the ferromagnetic transition, indicates a momentary 
decrease of the depletion rate. 

The free energy per electron and the many-electron energy contributions per electron 
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FIG. 7. (a) Total free energy per electron and (b) Hartree (solid), exchange (dashed) and correlation 
(dotted) energies per electron as a function of the held B. Vertical dashed lines indicate the 
transition fields in Fig. [5j 

are calculated dividing the corresponding magnitudes per unit length by the total electron 
density and plotted in Figs. [7](a) and (b), respectively. All energy contributions show weak 
kinks at Bp and Bl- The free energy per electron increases with B due to the increase 
in magnetic confinement, ffowever, the ffartree energy per electron (Fig. [7](b)) is reduced 
with B due to field-induced charge depopulation. At high magnetic fields, B > 18.5 T, the 
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Hartree energy changes sign because the free electron density is lower than the total density 
of static donors included in the simulation in the NW GaAs core. Note from Fig. [7^b) that 
the direct Hartree energy is one and two orders of magnitude larger than the exchange and 
correlation contributions, respectively, and therefore, it will rule many-electron effects in the 
system. 



FIG. 8. Non-interacting MSSs edges with respect to the InAs conduction band edge. Blue (red) 
dots are used for d-MSSs (d-MSSs) 

In order to asses the role of many-electron contributions, in Fig. [8] we show the MSSs 
calculated in a non-interacting model, i.e., vh — 0 and vxc = 0. The MSSs follow in this case 
a smooth evolution with B, which evidences the many-electron origin of the two transitions 
at Bp and B L in the SDFT calculation. We have also checked that such transitions persist 
when only vxc — 0 (data not shown here) as was expected from the weak effect of the 
exchange and correlation contributions in the present system (see Fig. [T)(b) ). 

Indeed, the transitions at Bp and Bl result from the balance between the two main 
energy contributions: the magnetic confinement, which increases the system energy with 
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B, and the direct Coulomb or Hartree energy, which is reduced with B due to the charge 
depletion decreasing in this way the system energy. Thus, the first transition at Bp, which 
produces an inversion in the slope of the |-MSSs, can be understood as a transition between 
a regime, B < Bp, in which the reduction in Hartree energy dominates over the magnetic 
confinement, to another regime, B > Bp, in which the magnetic confinement dominates. 
The key difference before and after Bp is the magnitude of the Hartree energy that is lost 
per depleted state, which is larger at B < B P . This is because when the system is not 
spin-polarized the Hartree energy also arises from the interactions between electrons with 
anti-parallel spin. The latter, which are absent in the ferromagnetic phase, are stronger than 
interactions between parallel spin electrons due to the lack of Fermi hole. 

The transition at B L , which produces an abrupt increase of the MSSs, is also interpreted 
with similar arguments, i.e., the Hartree energy lost per depleted state is lower at B > Bl- 
This is due to the larger localization of the electron density at B > Bl (cf. Figs. |5^e) and 
(f)), which entails a larger Fermi hole in the direct Coulomb interaction in this regime. 
Indeed, it has been proved that the conditional probability of finding an electron with a 
given spin when there is already another electron with the same spin nearby is lower when 
the former is localized. 35 

The spin-projected magnetoconductances G- j-, and the total magnetoconductance G = 
G^ + G± calculated from the SDFT modeling are shown in Fig. [9j^a) . Starting from low fields, 
the total magnetoconductance oscillates, due to oscillating MSSs crossing Ep, around an 
average value of 16 e 2 /h up to a magnetic field B ~ 6 T. As the held approaches Bp, a 
sudden step-like reduction of 4 magnetoconductance units is caused by the sudden depletion 
of the lowest set of f-MSSs (see Fig. [5](a)) induced by the ferromagnetic transition. 

At B > Bp the magnetoconductance shows an almost hat plateau that lasts up to B ~ 12 
T. This originates in the location of Ep in the symmetry induced energy gap between the 
second and third group of six f-MSSs (see Fig. [5] (a)). As Ep merges in the second group 
of J,-MSSs the magnetoconductance starts to oscillate again, while reducing in average at 
an increasing rate approaching Bl. At B > Bl Ep lies in the wide energy gap between the 
hrst and second Landau bands, hence G is constant. 

Finally, G(B) drops to zero when the hrst incipient Landau band crosses Ep and the 
conduction band gets completely depleted. In Fig. |9](b) we also plot the magnetoresistance 
1/G(B) to illustrate the kink observed at Bp, which corresponds to that observed in exper- 
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FIG. 9. (a) Total magnetoconductance G = (black) and spin resolved magnetoconduc¬ 

tances, G (blue) and G± (red), (b) Total magnetoresistance. Vertical dashed lines indicate the 
transition fields in Fig. [5j 

imental measures 36 37 of ferromagnetic transitions in flat quasi-ZD electron systems under 
in-plane magnetic fields. 


IV. SUMMARY AND CONCLUSIONS 

We performed a SDFT study of the electronic structure and magnetoconductance of 
hexagonal core-shell NWs pierced by an axial magnetic held. Critically, our modeling goes 
beyond often employed cylindrical and/or single-particle approximations to simulate radial 
hetero-structures, which neglect the strongly inhomogeneous, field-dependent distribution 
of the electron gas. 
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In the low field regime (B <■ 2T) we predict that AB magnetoconductance oscillations 
may disappear/resurface as a function of the gate-all-around voltage as a direct consequence 
of the presence of discrete symmetry induced energy gaps. Our calculations also allowed to 
critically analyze recent experiment^® 17 and justify the observation of AB oscillations in 
spite of the broken symmetry induced by the back-gate voltage. 

In the high magnetic field regime we found several field-induced transitions. First, the 
diamagnetic confinement induces a strong reshaping of the electron gas, which goes through a 
smooth transition from a low-field electron density distribution concentrated in the corners 
to a high-field distribution strongly localized in the facets of the radial heterojunction. 
Several experimental consequences of such reshaping are expected, for example in optical 
recombination experiments, due to the different localization of electrons and holes.® 

In addition, two abrupt transitions occur at discrete fields which are related to depletion 
of higher MSSs, either of the lowest anti-parallel spin MSSs, leading to spin polarization, or 
of the second incipient Landau-like band with parallel spin. The origin of these transitions 
lies in the increase of the effective Fermi hole occurring at each transition which affects the 
amount of Hartree energy that is lost per depleted state. As a consequence, such abrupt 
transitions are clearly marked in the calculated magnetoconductance by step-like behaviors. 
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